Lieb- Robinson Bounds and the 



Exponential Clustering Theorem 

Bruno Nachtergaele and Robert Sims 

Department of Mathematics 
University of California at Davis 
Davis CA 95616, USA 
Email: bxn@math.ucdavis.edu, rjsims@math.ucdavis.edu 

Version: August 12, 2005 

Abstract 

We give a Lieb-Robinson bound for the group velocity of a large class of discrete quantum 
systems which can be used to prove that a non-vanishing spectral gap implies exponential 
clustering in the ground state of such systems. 

1 Introduction 

One of the folk theorems in quantum lattice models claims the equivalence of the existence of a 
nonvanishing spectral gap and exponential decay of spatial correlations in the ground state. It has 
been known for some time that there are exceptions to one direction of this equivalence. There are 
models with a unique ground state with exponential decay of correlations but without a gap in the 
excitation spectrum above the ground state. For a simple example see Example 2 in 2, p 596]. In 
this paper we provide a rigorous proof of the other implication: a spectral gap implies exponential 
decay in the ground state. 

In relativistic quantum field theory it has been known for a long time that the existence of a mass 
(energy) gap implies exponential decay of spatial correlations. For example, in pQ Fredenhagen 
proves a general theorem applicable to arbitrary strictly local theories. It is the strict locality, 
i.e., the fact that space-like separated observables commute, not the relativistic invariance per se, 
which plays a crucial role in the proof of exponential decay. Non-relativistic models of statistical 
mechanics do not have strict locality, but there is a finite speed of propagation up to exponentially 
small corrections. This was first proven by Lieb and Robinson [li-t| . It is not a surprise that the 
Lieb-Robinson bound can replace the strict locality property. In particular, Wreszinski relied on 
it to prove a Goldstone Theorem in non-relativistic quantum statistical mechanics j^j. Precisely 
how to apply the Lieb-Robinson result to work around non-locality is not entirely obvious. Only 
recently, Hastings used it to derive exponential clustering for lattice models with a gap > an d to 
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obtain a generalization to higher dimensions of the Lieb-Schultz-Mattis theorem JU| in JTj. Our 
proof of exponential clustering is a rigorous version of Hastings' argument. 

In the proof by Lieb and Robinson the lattice structure played an essential role through the 
use of the Fourier transform. This was emphasized in a later version of their proof in |15| . Since 
lately there has been some interest in models without translation invariance or even without an 
underlying lattice structure such as spaces of fractal dimension [5] or the so-called complex 
networks |7||S], we revisited the Lieb- Robinson result and provide here a proof that works for quite 
arbitrary models defined on a set of vertices with a metric. 

2 Main Results 

We will consider quantum systems defined on a set of vertices V with a finite dimensional Hilbert 
space Tt x at each vertex x € V. At first we will assume that V is finite. For X C V, the Hilbert 
space associated with X is the tensor product Hx = te yH t , and the algebra of observables 
in X is denoted by Ax — B(7ix)- An interaction for such a system is a map $ from the set of 
subsets of V to A v such that $(X) e A x and $(A) = $(X)* for all X C V. The Hamiltonian is 
defined by 

H= £ $(A). 

XcV 

The dynamics of the model is the one-parameter group of automorphisms, {T t } tS R, defined by 

Tt (A) = e ltH Ae- ltH , A e Av- 

We will assume that V is equipped with a metric d. In the most common cases V is a graph, 
and the metric is given by the graph distance, d(x,y), which may be the length of the shortest 
path of edges connecting x and y in the graph. We will not, however, require an underlying graph 
structure for our results. In terms of d we define the diameter, D(X), of a subset X C V by 

D(X) = max{d(x, y) \ x, y € X}. 

Infinite systems can be introduced by considering a net of finite systems indexed by finite sets 
V partially ordered by inclusion. The C*-algebra of observables, A, is the norm completion of 
the union of the local observable algebras Av- We will assume that there is a uniform bound, N, 
on the dimension of the single site Hilbert spaces. An interaction <E> is defined as before but it is 
necessary to impose a boundedness condition in order for the finite-volume dynamics to converge 
to a strongly continuous one-parameter group of automorphisms on A. A standard reference for 
these infinite volume techniques is |15j . The strength of a given interaction will be measured by a 
norm, || ■ \\\, which for A > is defined by 

||$|| A := sup V \X\\\^(X)\\ N 2\x\ e *D(x) (1) 
xev xtx 

Here \X\ denotes the cardinality of the finite set X. For finite V, the supremum in this definition 
is of course a maximum. For infinite systems, finitcness of this norm is stronger than what is 
required for the existence of the dynamics, but it is this norm that appears in the Lieb-Robinson 
bound ^31- We will denote by B\ the set of all potentials for the system under consideration such 
that ||$||a < oo. 
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A Licb-Robinson bound is an estimate for the following quantities: 

C A>B (x,t) := [r t (A), B], (2) 

where x G V, t G R, A G A x , B G A. Due to the automorphism property of Tt, there is a symmetry 
in the norm of such commutators: || [r t (A), B]\\ = \\[A,T- t {B)]\\. 
It will be useful to consider 

C B (x,t):= sup — . (3) 

AeA* \\a\\ 

In a typical application we would have B G Ay, for some Y G V, and x G V\Y. Then, A and 
B commute and C B (x, 0) = 0. A Lieb- Robinson bound then aims to show that C B (x, t) is small 
for \t\ < T with T proportional to the distance between x and Y. 

Theorem 1 (Lieb- Robinson Bound) Fix X > 0, then for all $ G B\, x G V, t G K, and 

-B G -4, we /lave £/ie bound 

C B (x,t) < e 2 ^^C B (x,0)+ ]T e-^^^^e 2 ! 4 !!!*^-!)^^^) (4) 

It is straightforward to derive from Theorem ^ a bound for ||[rt(y4),B]|| for general local ob- 
servables A G Ax- One gets 

\\[T t (A),B]\\<N 3 W\\A\\'£c B (x,t) . 

For observables with overlapping supports, the Lieb-Robinson bound may not be useful in the 
sense that the trivial bound || [r t (A), B] \\ < 2\\A\\ \\B\\ may be better. The problem of estimating 
|| [rt (A), B) || for large t, in cases where it is expected to decay, is a separate issue that we do not 
address here. 

In the case of observables with widely separated supports the information in is essentially 
equivalent to 

C B ( Xl t) <Y,e 2ltlmx - Xd(x - y) CB(y,0) (5) 

For strictly local B G Ay, we can easily derive a more explicit bound by using that, C' B (y, 0) < 
2||£?||xy(2/), where xy is the characteristic function of the set Y. One obtains for all A G A x and 
B G Ay the bound 

\\[T t (A),B}\\<2\Y\\\A\\\\B\\e 2 W^- Xd ^ (6) 
If x (ji Y , a stronger inequality holds: 

\\[T t (A),B}\\ < 2\Y\ \\A\\\\B\\ (e 2 !*!!!*^ _ i) e -*d(x,Y) (7) 

It would be interesting to have an analogue of (JJJ for complex times. We may apply Theorem 
^with B replaced by th,(B), for any B which is analytic in a disk centered at of radius r > 0, 
e.g., a local B, and b < r. This way one can obtain, e.g., 

m^gjjCT < J^e 2 * l »- WM C Tft(fl) (»,0). (8) 

II II y£ y 
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The remaining open question is to obtain good bounds for C Ti6 (s)(y, 0). 

The next theorem provides an estimate on the spatial decay of correlations in states with a 
spectral gap. To state the gap condition precisely, we consider a representation of the system on 
a Hilbert space Tt. This means that there is a representation n : A —* B(TC), and a self-adjoint 
operator H on Ti such that 

7r(T t (A)) = e ttH ir(A)e- ltH , A e A. 

In the finite-system setting using a representation is, of course, merely for convenience and not 
essential. We assume that H > and that £1 £ Ti is a vector state such that Hil = 0. We say that 
the system has a spectral gap in the representation if there exists 6 > such that spec(i/)n(0, 6) = 
and in that case the spectral gap, 7, is defined by 

7 = sup{(5 > I spec(tf) n (0, 5) = 0}. 

Let Pq denote the orthogonal projection onto ker H . From now on, we will work in this represen- 
tation and simply write A instead of ir(A). 

We will derive a bound for ground state correlations of the form 

(n,A nb (B)Q) (9) 

where b > and A and B are local observables. The case b = is the standard (equal-time) 
correlation function. 

Theorem 2 (Exponential Clustering) There exists fj, > such that for all A 6 A x . B 6 Ay, 

x £ Y, for which PqBCI = PqB*Q — 0, and b sufficiently small, there is a constant c(A, B) such 
that 

\{Cl,AT ib {B)Cl}\ < c(A, B)e~' ld{x ' Y) ( 1+ ^<*,^) . (10) 

One can take 

c(A,B) = \\A\\\\B\\ (1+^1+ 1 ], (11) 

and 

^=4M^- (12) 

The bound is valid for < 76 < 2fid(x, Y). 

Usually, the distance between x and Y, will be bounded below by some a > 0, so c(A, B) in 
(lll|l can be replaced by a constant which depends only on the norms of A and B and the size of 
their supports, but not the distance between them. 

Note that in the case of a non-degenerate ground state, the condition on B is equivalent to 
(57, Bil) — 0. In that case, the theorem with 6 = becomes 

I (fl, ABil) - (fi, Aft) (O, SO) I < c(A, B) e -» d{x ^ . (13) 

For b > large, there is a trivial bound 

\{Q,Aru,{B)Sl)\ < \\A\\ \\B\\ e -^ b . (14) 

For small b > 0, the estimate <|10|) can be viewed as a perturbation of (|13l) . Often, the important 
observation is that the decay estimate l|10|) is uniform in the imaginary time ib, for b in some 
interval whose length, however, depends on d(x,Y). 
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3 Proofs 

3.1 Lieb-Robinson bound. Proof of Theorem [T] 

Our argument closely follows the proof of [151 Proposition 6.2.9], but we avoid the use of the 
Fourier transform in order to obtain a generalization to arbitrary sets of vertices with a metric. 

Let A £ A x and B £ A, and consider the quantities Ca,b{%, t) and C B (x, t), defined in J2J and 
©, respectively. 

From the fundamental theorem of calculus, one has that 



C AtB (x,t) = C AtB (x,0) + i I [t s ([$(X),A]), B]ds. 



Several applications of the triangle inequality to (|15J) yields 

"1*1 [r,([<I»(.Vl..i;i. J fJ] 
i sup - 

1 r, A r- A 

X3i' 



r\t\ 

C B (x,t) < C B (x,0) + V / sup 



ds. 



Now, for any finite X d V with x G X, one may write 



\x\ \x\ 

[*po, a] = y, E C U ■>• s ) II' a, :•• 

i=l j Xi = l,...,N*. »=1 



(15) 



(16) 



(17) 



where each Xi € X and e(j Xi ) is a matrix unit for the algebra A Xi ■ For each term in l|17|l. expanding 
the commutator of the product yields 



1*1 



Ur s (e(j Xi )),B 



vex 



(18) 



Combining the linearity of r s , the basic estimate 1)18[L and the fact that the coefficients in the 
above expansion H17|l satisfy 



we arrive at the inequality 



C(j xl ,...J xm )\ < 2\\<Z>(X)\\\\A\\, 

, r \t\ % 

C B (x,t) < C B (x,0) + 2 J2 ll$P0ll^ 2|x| / £C B (y, S )d S . 

X3x ^° »'CY 



Motivated by the expression above, we define the quantity 



and rewrite (|20|l as 



e(x,y) := £ ||$(X)|| A^l 



CbCi,*) < C B {x,0) + 2 V e(x,y)C B (y,s)ds 



(19) 



(20) 



(21) 



(22) 
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Iteration of l|22[l yields 

C B (x,t) < C B {x,0) + 2\t\Y e(x,y)C B (y,0) 



y ev 



Y < x >y) E <y,y')c B (y',o) 

yeV y'eV 



(2|t|) 



3! 

y ev y'ev y"ev 



Recall that we assumed that the interaction satisfies a bound of the form 

||* ||a = sup V \X\\\<£{X)\\N 2 \ x \e XD W < oc, (23) 
^ v itx 

for some A > 0. We wish to prove an exponential bound on the quantity C B (x,t). To make this 
bound explicit, we set e x (x,y) := e xd ( x,y ^ e(x,y), and observe that 

= E E MX)\\N 2 ^e Xd ^ (24) 

yeV yeV X3x,y 

< J2 E UWII^e^W - ||$ ||a. 

X3x yeX 

Now, returning to the iterated version of 112211 above, we may resum the upper bound, which is 
allowable as all the terms are non-negative. Using the triangle inequality often, one derives 

C B (x,t) < e- Xd ^f(x,y)C B (y,0), (25) 

yeh 



where 



f(x,y) = 6 x , y + 2\t\e x (x,y) + Y e x (x,y')e x (y',y) 



2 

y'ev 



(2|i|) 



Y E e x (x,y")e x (y",y')e x (y',y) 



3! 

y "ev y'ev 

Moreover, summing over all but the first term appearing in / and using 124|) . we see that 

f(x,y) < S x , y + ( e 2 l*lll*H>-l). (26) 
Therefore, we have proven the estimate 

C B (x,t) < e 2 1*1 11*11* C B (z,0) + Y e~ Xd ^ (e 2 1*1 " * «> - l) C B (y,0) , (27) 

y€V: 

yjtx 

as claimed. □ 
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3.2 Decay of correlations. Proof of Theorem [2] 

It is well-known, see e.g. ^1^1 > that for interactions <& S B\ and B a local observable the function 
1 1 ► T t (B) is analytic in a disk centered at the origin. 

Part of our argument requires analyticity for arbitrarily large \t\. Therefore, we consider entire 
analytic elements B(a), a > 0, (see Proposition 2.5.22]), with the property that ||B(a)|| < ||£?|| 
and lim a ^o \\B(a) — B\\ =0. In this case, the function 

/„(*) := (n,AT z (B(a))Cl) (28) 

is entire. Observe that the quantity of interest, i.e. that which appears on the left hand side of 
(I10|) . corresponds to fo(ib) where /o is as in (|28(l above with B(a) = B and b > 0. The case b = 
follows by a limiting argument. 

Take T > 2b and denote by Tt the contour in the complex plane which joins — T to T along 
the real axis and T to — T via the semi-circle Te l6 for 9 6 [0,7r]. By analyticity, we have that 

U(ib) = [ ^rdz. (29) 
tow JT T z — ib 

Using the existence of a gap 7 > and the assumption PqBVL = 0, one easily finds that, for 
0€[O,7r], 

l/a(7V S )| = 

From this, we conclude that 

|/ (i6)| = lim|/ a («OI < 

a— >0 

where we have used that ||B(a) — B\\ —> as a — + 0. Since the second term vanishes in the limit 
T —>■ 00, we have shown that 

\(ft, AT ib (B)Q)\ < limsup 

and thus, the proof has been reduced to estimating this integral over the real axis. 

Our estimate follows by splitting the integrand into three terms. Let a > and write 

fo(t)e a " 2 = f (t)e- at ' 2 + f (t) (e ab2 - e""* 2 ) . (33) 

The first term above may be written so that the commuator explicitly appears, i.e., 

fo(t)e- at2 = (n,r t (B)An)e- at2 + (S!,[4,T f (B)]!l)r« i! . (34) 

Using (|33|l and (|34|) . we arrive at a bound of the form: 

< I + II + III , (35) 



n,Ae lTeZ<>H B(a)n) < \\A\\ \\B\\e 



- T 7 sin ( 



(30) 



1 

2~7ri 



— T 



fo(t) 
t - ib^ 



dt 



\M\\B\ 



-Tjsin 



(31) 



1 

2~7ri 



t - ib" 



dt 



(32) 



1 

2~7ri 



— T 



hit) 

t - ib" 



dt 
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where I corresponds to the intergral containing the first term of (|34|) , 77 contains the second term 
of (|34|l . and III the second term of (13311 . Each of these terms will be bounded separately in the 
limit T — > oo. 



The first term: The first term appearing in the bound on the right hand side of i|35[l is 

- f 

2ni J _ T 

Using the spectral theorem, we may write 



i_ f 1 (n,T t (B)An)e- at - 

t-ib 



dt 



(Q,T t (B)AQ) = (e ltH B*n,Ail) = / e - itE d(P E B*n, Aty 7 



(36) 



(37) 



where Pe is the projection valued spectral measure corresponding to H, and we have used that 
PoB*tt = 0. One may now write 



(n,Tt(B)An)e- at ^ dt _ r°° f T e~ UE e~ ai ' : 



t - ib 



t-ib 



dtd(p E B*n,An). 



/-y J —T 

The inner intergral can be estimated using Lemma ^ proven below, from which we have 

hm sup I < - — — — - e 4 ° , 

independent of b. 



(38) 



(39) 



The second term: The second term is the one for which we shall apply the Lieb-Robinson estimate 
in the form of 0: 



II := 



j_ r T (n,[A,T t (B)]n)e- at ~ 

2iri I _rp t — ib 



dt 



arifie-^,, (40 

" 2W-00 \t\ 



Here it is crucial that the "time", ib, is purely imaginary. We will now break the integral up into 
two regions: {t € M : \t\ < s} where the Lieb-Robinson bound is useful and {t € M : \t\ > s} in 
which we will use a norm estimate: 



1 

2^ 



[AMB)\ 



dt 



^ 2\Y\ \\A\\ \\B\\ cMmw _ M{xX) + \\A\\\\B\\ e _ as ^ 
7T s J ix a 



(41) 



The third term: The final term is 



III := 



2ni 



(il,A Tt (B)Sl) (e c 



-T 



t-ib 



-dt 



(42) 
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.2 

which represents the penalty we incur for introducing the gaussian factor e~ and the scaling 
e ab . This term can be made small by a judicious choice of parameters. Using the spectral theorem 
again, we have 



III = 



1 

2ttI 



t e l 

-T 



t~ib 



dtd{A*n,P E BVL) 



(43) 



By adding and substracting the appropriate limiting quantitites, which are shown to exist in 
Lemma ^ the integrand above can be rewritten as 



af -Eb 

e e 



1 



e 4 ° 



dw + e ab ~ (F ib<T (E) - e- Eb ) + i?i + i? : 



(44) 



The quantities Fh,,t,Ri, and i?2, are defined in lf^T|l . (|ST|) . and respectively. Using 1(551) for 
complex t, one may estimate the first term in (|44|) as 



2 A /7ra 



e - bw e —dw < -e-~° 



(45) 



for 2ab < 7 and E > 7. Using dominated convergence, one sees that the other terms vanish in the 
limit T — > 00. This implies the following bound: 



limsup/n < ~\\A\\ \\B\\e~ 

T— »oo 2 



The proof is completed by choosing a = 3- and s so that 

s(2\\4>\\ x + l) = Xd(x,Y). 



We have that 



\(n,An b (B)n)\ < \\A\\ \\B\\ 1 



2 in 



where 



7T A / 7 r//d(a;,y) / 
A 7 



" 4||#||x + 7' 

This proves Theorem [5] 

In the proof above we used the following lemma, which is a variation of Lemma 3.1 in |17j 

Lemma 1 Let E E R, a > 0, and z 6 C + := {z e C : Im[z] > 0}. One /ias ttai 



(46) 
(47) 

(48) 

(49) 
□ 



lim 

T^oo 27TI 



_ T i — z 

Moreover, the convergence is uniform in z G C + . 



e™* e-^r- dw. 



(50) 
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Proof. One may prove this lemma by observing some simple estimates on the following function: 



2m J_ 



t - z 



dt. 



(51) 



For w > 0, integration over the rectangular contour Tt joining — T —> T ^ T+iT — > —T+iT — > — T 
yields 



1 



dz' = 



2ni J Vt z' — z 
for any < 2\z\ < T. From this one may conclude that 



\F Z}T (w) - e iwz \ < - 



1 

wT 



1-e 



-wT\ 



-wT 



If w < 0, then closing a similar rectangular contour in the lower half plane yields 

2 



\F z ,tH\ < 



\w\T ( 



1 - e 



\w\T 



\w\T 



(52) 



(53) 



(54) 



Now, as the Fourier Transform of a gaussian is another gaussian, one easily sees that for any 
E e K and a > 0, 

' ~^dw. (55) 



2v7ra 



e lwt e 



We may now write the prelimit quantity on the left hand side of (|50|l as the desired limit and two 
remainder terms, i.e., 



1 

2ttI 



T e iEt e -at 2 



dt 



1 



2yhf^ 

1 

2~Jna 



e 4 Q F Z: T(w)dw 



e ™* e -iH_a_ rfiy + Ri + R2 



(56) 



Using the estimates (1531) and (|53|l. one may apply dominated convergence to conclude that both 
the remainders 

1 /-00 2 

Hi := tt^= / e-*^ 1 - (F ZtT (w) - e lwz ) dw (57) 



2,h^ 



and 



vanish in the limit T — > oo. 



1 



_ (vi-E)' 1 

e 4c« F Zi T(W)dw 



(58) 
□ 



4 Applications and Generalizations 

For many applications it is of interest to consider the thermodynamic limit. Our clustering theorem 
allows for two approaches to the thermodynamic limit, each with its own merits. In the first 
approach, one applies the theorem to finite systems and obtains estimates that are uniform in the 
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size of the systems. This may require a careful choice of boundary conditions. The bounds then 
carry over automatically to any state which is a thermodynamic limit of finite volume states for 
which the estimates are obtained. This is the most straightforward way to proceed in cases where 
the finite volume ground states are unique, weak* convergent and with a uniform lower bound on 
the spectral gap. 

In the second approach one focuses directly on the infinite system. This requires a proof of 
a spectral gap for the infinite system, but may avoid difficulties associated with boundary states 
which may obscure the spectral gap for finite volumes. 

Another method for dealing with boundary states is to generalize Theorem [5] to systems with 
quasi-degenerate ground states. As long as the quasi-degenerate states remain separated from 
the excitation spectrum by a non-vanishing gap, this is a straightforward generalization, but the 
condition PqBSI = 0, where now Pq is the spectral projection on all states below the gap, may be 
a bit tricky to verify. 

Another rather obvious generalization of our results is to systems of fermions on (V,d). For 
even interactions (only products of an even number of creation and annihilation operators occur), 
and even observables, all results carry over without change. If A and B are both odd observables, 
one gets a Lieb-Robinson bound for the anticommutator instead of the commutator. 
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